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ABSTRACT 

In 2009, Cygnus X-3 (Cyg X-3) became the first microquasar to be detected in the GeV 7-ray regime, via 
the satellites Fermi and AGILE. The addition of this new band to the observational toolbox holds promise for 
building a more detailed understanding of the relativistic jets of this and other systems. We present a rich 
dataset of radio, hard and soft X-ray, and 7-ray observations of Cyg X-3 made during a flaring episode in 
2010 May. We detect a ^3-d softening and recovery of the X-ray emission, followed almost immediately by a 
~1-Jy radio flare at 15 GHz, followed by a 4.3er 7-ray flare (E > 100 MeV) ~1.5 d later. The radio sampling 
is sparse, but we use archival data to argue that it is unlikely the 7-ray flare was followed by any significant 
unobserved radio flares. In this case, the sequencing of the observed events is difficult to explain in a model 
in which the 7-ray emission is due to inverse Compton scattering of the companion star's radiation field. Our 
observations suggest that other mechanisms may also be responsible for 7-ray emission from Cyg X-3. 
Subject headings: black hole physics — X-rays: binaries — X-rays: individual (Cygnus X-3) 



1. INTRODUCTION 

The X-ray binary (XR B) system Cyg n us X- 3 (hereafter 
Cyg X-3), discovered by iGiacconi et al.1 (119671). is notable 
for it s Wolf-Rayet (WR) comp anion (Ivan Kerkwiik et alJ 
1996), short (4.8 h) orbital period dBonnet-Bidaud & C hardin 
1988), and dramatic radio variability (Waltman et al. 1994). 
The relativistic jets that it produces (Geldzahler et al.lll983l) 
classify the source as a "microquasar" system. The dis- 
tance to Cyg X-3 is ~7 or ~9 kpc, depending on the 
method used dPredehl et al.1 120051 iLing et al.1 120091) . It is 
also uncertain whether the compact object is a black hole 
(e.g. [C herepashchuk & Moffat] 119941) or neutron star (e.g. 
iStark & Saiall2003l) . 

Cyg X-3 became more notable in 2009 when it be- 
came the first microquasar to be detected in the G eV 7-ray 
regim e, via the satellite s Fermi dAbdo et al J 120091) and AG- 
ILE dTavani et alJ 12009b . A 4.8-h modulation in the Fermi 
data makes the detection definitive. Emission of 7-rays has 
been detected by both observatories during every "high/soft" 

'Department of Astronomy, 601 Campbell Hall # 3411, 
University of California, Berkeley, CA 94720-3411, USA; 
pwilliams@ astro .berkeley.edu 

2 Space Sciences Laboratory, 7 Gauss Way, University of California, 
Berkeley, CA 94720-7450, USA 

3 The University of Cambridge, Milliard Radio Astronomy Observatory, 
Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CB3 0HE, UK 

4 CRESST and NASA-Goddard Space Flight Center, Astrophysics 
Science Division, Astroparticle Physics Laboratory, Greenbelt, MD 
20771, USA 

5 CSST, University of Maryland Baltimore County, 1000 Hilltop Circle, 
Baltimore, MD 21250, USA 

6 AIM-Astrophysique Instrumentation Modelisation (UMR 7158 
CEA/CNRS/Universite Paris 7 Denis Diderot), CEA Saclay, 
DSM/IRFU/Service d'Astrophysique, Bat. 709, L'Orme des Merisiers, 
FR-91 191 Gif-sur-Yvette Cedex, France 

7 Dr. Karl Remeis-Sternwarte and Erlangen Center for Astroparticle 
Physics, Universitat Erlangen-Numberg, Sternwartstr. 7 , 96049 Bamberg, 
Germany 

8 Departament d'Astronomia i Meteorologia (DAM) and Institut de 
Ciences del Cosmos (ICC), Universitat de Barcelona, Martf i Franques 1, 
08028 Barcelona, Spain 

9 Special Astrophysical Observatory RAS, Karachaevo-Cherkassian 
res, Nizhnij Arkhyz 36916, Russia 



X-ray state to occur since they began science operations 
(Abdo "etani2009l iTavani et al.l[2009t iCorbel & Havsll2010l 
®. The high/soft state is also associated with radio flaring 
and is closely tied to the launching of relativistic jets (e.g. 
iFender et al.ll2004t iKolionen et al.ll2010l) . 

Various physical processes are expected to be involved 
in the emission and propagation of high-energy photons in 
Cyg X-3. Inverse Compton (IC) upscattering of the intense 
(~10 5 erg cm" 3 ) radiation field of the WR companion by 
high-energy (7 > 10 3 ) electrons in the re lativistic jet is the 
most efficient source of GeV emission dAbdo et al] |2009; 
Dubus et al. 2010). If the density of relativistic protons in the 
jet is sufficiently high, inelastic collisions between these and 
protons in the stellar wind will produce neutral pions that de- 
cay in to 7-rays (among other products) as well (R omero et al] 
2003). In this case, there will be further 7-ray emission due 
to secondary leptons created during the hadronic interactions 
dOrellana et al.ll2007b . The typical hadronic content of mi- 
croquasar jets is unknown: large numbers of hadrons impose 
significant energetic constraints, but there is dire ct evidence 
of their presence in at least one system, SS 433 dKotani et al] 
[1991 [19961: iMigliari et al.ll2002l) . 

A 7-ray flare due to IC upscattering by jet electrons will 
be followed by a radio flare a few days later as those electrons 
cool, rarify, and become transparent to their radio synchrotron 
emission. This phenomenon has been observed rep eatedly but 
not invariably dAbdo et al]|2009t ITavani et al.ll2009l) . If the 7- 
ray emission is primarily hadronic, the result will be substan- 
tially the same, as secondary p rocesses will creat e a similar 
population of energetic leptons. lAbdo et ail {2009) determine 
an approximate 5 ± 7 d lag between radio and 7-ray emis- 
sion. Inspe c tion o f their Fig. 2 and the data presented by 
Tavani et al. (2009) suggests that this uncertainty is primar- 
ily systematic rather than statistical; i.e., that the radio/7-ray 
lag varies significantly from event to event. Blazar observa- 
tions often show lag variability, both between sources and in 
the same source over time, with theoretical supp ort for sign 
chang es in the lag between 7-ray and other bands (All er et ail 

Hoii. 

In 2010 May, multiple monitoring programs detected 
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that Cyg X-3 was becoming more active in 7-ra y s and 
was entering the high/soft state dBulgarelli et al.l 
Corbel & Fermi Large Area Telescope C ollaboration 
Kawai et alJ 120101: ICorbel & Haysi 120101) . We describe 
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multiband observations of Cyg X-3 made during this episode 
and present our results (§|3}. Our data show a 7-ray flare 
apparently lagging a radio flare, a sequencing opposite of 
that expected in the typical interpretation. We discuss and 
interpret this result (§@) and, finally, present our conclusions 

2. OBSERVATIONS & DATA ANALYSIS 

We observed Cyg X-3 during its 2010 May flaring episode 
with a variety of instruments. Below, we describe the obser- 
vations and their analysis. Many dates and times we mention 
are near MJD 55,340. For brevity, we express these in an 
offset MJD, defining OJD = MJD -55, 300. Thus OJD is 
JD 2,455,300.5 or 2010 Apr 14.0 UT, and the 7-ray flare was 
detected on OJD 43 = 2010 May 27 UT. 

2.1. Radio 

The Arcminute Microkelvin Imager (AMI) arrays (Cam- 
bridge, UK) are two aperture synthesis telesco pes mainly 
used to study the cosmic microwave background (Zwa rt et all 
2008). The observations described herein were made with 
the Large Array (AMI-LA), the reconfigured and reequipped 
Ryle Telescope, consisting of eight 13-m antennas with a 
maximum baseline of about 120 m, observing in the band 
12-18 GHz. The angular resolution is typically 25". Mon- 
it oring of small-diameter s ources is undertaken as described 
in Poolev & Fender! (11997) ; observations are interleaved with 
those of a phase-reference calibrator, and after calibration the 
data for individual baselines are vector-averaged. The in- 
phase component then provides an unbiased estimate of the 
target source's flux density. The amplitude scale was cali- 
brated by (at least) daily observations of 3C 48 and 3C 286, 
both of which are believed to be very nearly constant on long 
timescales. 

Should there be emission on a scale resolved by some 
of the baselines, the flux density estimate would be incor- 
rect. There is an extended region of low- brightness emission 
around Cyg X-3 (Sanchez-Sutil et al. 2008) which is detected 
on the shortest baselines at extreme hour angles with the AMI- 
LA, and is sufficiently bright to need separate treatment when 
the central source in Cyg X-3 is "quenched" (flux density < 
10 mJy at 15 GHz). No such correction was required in the 
observations described here. 

Observations of Cyg X-3 were also made with the RATAN- 
600 telescope (Nizhnij Arkhyz, Russia) as part of an on- 
going microquasar monitoring campaign (Trushkin 2000; 
Trush kin et alJ [2006). The o bservations and data analysis 
were performed as described in Trushk in" et al.l(t200 6). In this 
work we show results at 1 1 .2 GHz; near-simultaneous obser- 
vations were made at 4.8 and 7.7 GHz as well. 

Cyg X-3 was observed b y the Allen Telesco pe Array (ATA; 
Northern California, USA; IWelch et al.ll2009b six times in the 
period OJD as part of a larger transient search. Con- 

tinuum images at 3.09 GHz were made after subtraction of 
model of the static sky. Source fluxes were determined in 
— 10-minute segments by fitting point source models in the 
image domain. 

2.2. X-Ray 



Our study includes hard X-ray observations made with the 
International Ga mma-Ray Astrophys ics Laboratory (INTE- 
GRAL) satellite dWinkleret all 120031) . We used the INTE- 
GRAL Soft Gamma-Ray Imager (ISGRI: iLebrun et alj|2003h 
instrument and included two observations that were part of 
a program to observe Cyg X-l (PI: Wilms). Although the 
primary target of the observations was Cyg X-l, the ISGRI 
FOV is large enough to include the entire Cygnus region, and 
our target was in the FOV throughout the observations. The 
first observation occurred during revolution 929, starting at 
OJD 40.625 and ending at OJD 41.940. The second observa- 
tion occurred during revolution 938, starting at OJD 67.469 
and ending at OJD 68.8 05. A preliminary re port of the first 
observation was given in lTomsick et al.l (120101) . 

We reduced the ISGRI data using the Off-line Scien- 
tific Analysis (OSA v9.0) software package. The program 
ibis_science_analysis is the primary tool for ex- 
tracting the standard data products. We produced and in- 
spected the ISGRI image in the 20-40 keV band, and then 
made energy spectra and light curves for all of the bright 
sources in the field. After correcting for instrumental dead- 
time, the total exposure times for Cyg X-3 from the first and 
second observations are 62,420 s and 66,540 s, respectively. 
We produced 20^-0 keV light curves with a time resolution 
of 1 ks. 

We also obtained one-day average quick-look measure- 
ments from the All-Sk y Monitor on the R ossi X-Ray Timing 
Explorer (/^TfiyASM ^Levine et aljl996l) and transient mon- 
itor results from the Burst Alert Telescope on the Swift 7-ray 
burst mission (Swift/BXT; Barthel mv et al.ll2005l) . These re- 
sults are provided publicly by the RXTEI 'ASM and Swift/BKT 
teams, respectively. 

2.3. "/-Ray 

Fermi Science Tools 9.17 and H EASoft 6.9 were us ed to 
reduce and analyze all Fermi-L AT dAtwood et alj|2009t) ob- 
servations within 10° of Cyg X-3 that took place in the range 
OJD -59-150 (2010 Feb 14 - 2010 Sep 1 1 UT), providing a 
baseline of -100 d before and after the flare of OJD 43^44. 
Cyg X-3 is within the LAT FOV for — 15 ks per day. Given 
this and the typical flux of Cyg X-3 in the LAT band, we were 
compelled to use 1-day time bins in the following analysis. 

Thirty arcminutes away from Cyg X-3 is a compara- 
tively bright pulsar, PSR J2032+4127. Using gtephem 
and gtpphase, we ad ded the phase fr om the most re- 
cent ephemeris available (Abdo et al. 2009) to the events file. 
Then, we extracted events correspond ing to the off-puls e 
phases (0-0.12, 0.2-0.6, and 0.72-1; ICamilo et all [2009h . 
This removed 20% of the live time, so the exposure time of 
the off -pulse events file was corrected accordingly. 

With gtselect and gtmktime, good time intervals 
from the off-pulse events file were selected in the 0. 1-10 GeV 
energy range. At high energies (>10 GeV), Cyg X-3 is ex- 
pected to have negligible emission. The events class was 
set to "3" which selects only high-quality diffuse-background 
photons. To minimize background albedo photons from the 
Earth's limb, zenith and rocking angles were restricted to less 
than 105° and 52°, respectively. 

Exposure maps were generated by gtexpcube and 
gtexpmap while gtbin created photon counts maps in the 
region of interest. Emission from Cyg X-3 was not apparent 
on the counts maps, which is unsurprising given that Cyg X-3 
is a faint source whose peak emission is at the low end of the 
LAT band where the point spread function is 5-10°. Compli- 
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eating matters is the high level of diffuse background emis- 
sion. 

Following the procedure described in Ab do et aTl (|2009), 
unbinned likelihood analysis was performed for the daily 
bins with gtlike consideri ng photons inside a 7° radius 
of Cyg X-3 from all 1FGL (lAbdo et all 120101) sources up 
to 5° away, and all bright (detection significance > la and 
flux [>100 MeV] > 5 x 10" 8 ph cm -2 s _I ) 1FGL sources up 
to 20° away. The instrument response function was "Pass 
6 v3" (P6_V3_DIFFUSE), and the convergence relied on 
the NEWMINUIT method. Spectral models and input param- 
eters from the 1FGL catalog were used except for pulsars 
1FGL J2021. 0+3651 and 1FGL J2021. 5+4026 which had ex- 
ponential cutoffs in addition to their power laws. Besides 
Cyg X-3, whose photon index and normalization were free, all 
sources had spectral parameters fixed to the 1FGL values. We 
included models for the galactic (gll_iem_v02) and ex- 
tragalactic (isotropic_iem_v02) diffuse emission. The 
normalizations were left free to vary, as was the photon index 
of the galactic component. 

3. RESULTS 

Long-term lightcurves of Cyg X-3 in soft and hard X-ray 
bands are presented in the top panel of Figure Q] The system 
was in a soft X-ray state during OJD ~20-55, with an episode 
of particular softness occurring during OJD 39—41. 

The Fermi source flux integrated over 0.1-10 GeV is pre- 
sented in the lower panel of Figure Q] In bins where the test 
statistic (TS ~ -21nL, where L is the ratio of the likelihood 
of models without and with the source) is less than one, we 
compute and plot a la upper limit. A peak is seen starting on 
OJD 43. The TS value in this bin is 18.4, which translates ap- 
proximately to 4.3<r, and the flux is (4 ± 1) x 10~ 6 ph cm" 2 s" 1 . 
We are unable to achieve a more precise timing of the flare. 
Shifting the binning by half a day results in lower-significance 
and lower-flux detections in the range OJD 42.5^4.5. Anal- 
ysis using higher-cadence binning does not yield significant 
detections of the peak. 

Radio measurements from just before the 7-ray flare are 
presented in the top panel of Figure |2] A ~1 Jy radio flare oc- 
curred on OJD 42.25. There is evidence of an earlier flare 
on OJD 41.11, with a maximum observed flux density of 
366 mJy, but only the rising portion of the flare is observed, 
so its properties are poorly-constrained. As is shown in Fig- 
ure [3] the 15 GHz flux density of Cyg X-3 is depressed in the 
seven days following the 7-ray flare, with a mean flux den- 
sity of 48 mJy and an 80 th -percentile observed flux density of 
66 mJy. (That is, 80% of the measurements during this pe- 
riod are < 66 mJy.) Starting on OJD 51, the typical radio flux 
density increases by a factor of ~2 to ~100 mJy. 

The bottom panel of Figure [2] shows high-cadence hard 
X-ray measurements from before the 7-ray flare, made as 
Cyg X-3 exited the softest portion of its soft state. They show 
the recovery of the hard X-ray flux, subject to orbital mod- 
ulation. The modulation makes it difficult to determine the 
significance of the drop in flux seen in the last three mea- 
surements; it may be evidence of a significant reduction in 
hard X-ray emission precursing the radio flare. Combining 
the Swi/f/BAT and INTEGRAL measurements, the hard X-ray 
flux has recovered by OJD ^42.0, about 5 h before the ob- 
served radio flare. The 7-ray flare then occurred ~l-2 d later. 

4. DISCUSSION 



If our observations are taken at face value — that is, one as- 
sumes no significant radio activity during coverage gaps and 
that the flaring stems from one ejection event — the sequenc- 
ing of the radio and 7-ray emission is inconsistent with the 
companion IC model discussed in <JT] Emission in 7-rays 
following a radio flare could be interpreted as a reenergiza- 
tion o f the relativistic jet electrons by a shock (Dub us et al. | 
201C ). With a typical jet speed of ~0.5 c dTudose et al.l 
2007) and a delay of ~1 d, the reenergization would occur 
^100 AU from the system. This is much larger than the or- 
bital separation but three orders of magnitude smaller than 
the distance at which the jet is expected to interact signif- 
icantly with the ISM (~1 pc). Such reenergization could 
be due to a co llision between the jet and a dense clump of 
the WR wind (Araudo et al. 2009). In this scenario, the ab- 
sence of 7-ray emission at the time of ejection can be ex - 
plained by absorption within the system (ICerutti et al.112.01 ll) . 
alt hough this model m ust b e reconciled with th e observations 
of lTavani etail d2009l) and lAbdo et all d2009l The reener- 
gization would lead to additional radio emission as the 7-ray- 
emitting electrons cool, although the ambient magnetic field, 
and hence synchrotron luminosity, will be much weaker than 
that found close to the system. 

If the 7-ray flare is interpreted as the result of a dis- 
crete ejection event, the lack of a notable subsequent radio 
flare could be explained by the ejected material being largely 
hadronic. While hadronic 7-ray emission is not as efficient 
as IC upscatteri ng, it does yield a h igher ratio of 7-ray to ra- 
dio luminosity (Rom ero & Vilall2008l) . The secondary leptons 
due to hadronic interactions would, however, radiate, and de- 
tailed studies typically find that their bolom etric luminosity 
is com parable to that of the primaries (e.g. IVila & Romerol 
l2010h . 

Should our data be taken at fac e value? App lication of 
a simple synchrotron-cloud model (Ivan der L aan 1966J) sup- 
ports the intuition that radio flares with sizes comparable to 
that of the largest observed could have occurred without de- 
tection, with the modeled lifetimes being ~0.2 d. While this 
m odel has had ambiguous success when applied to Cyg X- 
3 (Fen der et al.l 11997). its simplicity is advantageous for our 
sparse data, and it has been successfully applied to observa- 
tions of other systems (e.g. iFender et al.l 1 19971 IWilms et al.1 
120071; iPrat et alJl20Tob . To obtain a better understanding of 
the radio behavior of Cyg X-3 we compared the 2010 May 
data to 15 years of archived AMI-LA / Ryle Telescope obser- 
vations. The ~1 Jy radio flare is unusual: no flux densities 
> 400 mJy are detected in a 500 d span around the observa- 
tions we describe. 

In light of the predictions of the companion IC model, we 
consider particularly the radio observations after the observed 
7-ray flare. The duty cycle of the observations is only 6.8% 
in the seven-day post-flare quiescent interval. We identified 
5 1 epochs similar to this in the archives, each lasting at least 
7 days and having an 80 th -percentile observed flux less than 
70 mJy (cf. ®. These archival measurements have a duty 
cycle comparable to that of the 2010 May data (4.9%), but a 
much larger time on-source (55.5 days). Because this time is 
large compared to the total duration of the 2010 May post- 
7-flare quiescence, a search for flares in the archival dataset 
can constrain the likelihood of there having been an unob- 
served flare in the 2010 May dataset, assuming no long-term 
evolution in the behavior of the quiescent state and stochastic 
flaring. 

The epochs were identified by exhaustively searching for 
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Figure 1. Long-term lightcurve of Cyg X-3 in X- and 7-rays. Shaded region: timespan depicted in Figure[2] Vertical black line: time of maximum observed 
15 GHz flux density. On some days, our derived Fermi upper limits are smaller than the canonical background rate, 0.3 X 10~ 6 ph cirT 2 s ; on those days, we 
use the background rate as the upper limit. 
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Figure 2. High-cadence lightcurve of Cyg X-3 in radio and hard X-rays in the time immediately before the 7-ray peak (= OJD 43 -44). The X - rays a re subject 
to orbital modulation with a 4.8-h periodicity. Vertical blue markers: times of X-ray minima according to parabolic ephemeris of Singh et al. (2002). Shaded 
region : partial timespan of Fermi peak. 



seven-day segments meeting the aforementioned criteria, then 
lengthening these segments as much as possible without vio- 
lating the 80 th -percentile constraint. Segments separated by 
less than three days were merged, in a few cases shortening 
the segment somewhat to preserve the statistical constraint. 
This method of construction does not bias against epochs con- 
taining rapid flares, no matter their size. Qualitatively, the 
lightcurve of the typical epoch starts high, drops to very low 
flux densities, and then becomes high again, possibly with 



rapid flares in the middle. 

The largest flux density seen in the selected archival data 
was 524 mJy, in the context of a single rapidly-evolving (rise- 
time ^0.1 d) flare that would not have necessarily been de- 
tected in the 2010 May observing. (Here we exclude slow, 
large, epoch-terminating flares that would have been eas- 
ily seen in the 2010 May data.) Approximating the rate of 
such flares as one per 55.5 days of observing, we derive an 
11% chance that such a flare occurred in the 7 days after the 
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Figure 3. Radio lightcurve of Cyg X-3 during and after the 7-ray flare. The three panels proceed in chronological order from top to bottom. Red arrow: time 
of maximum observed 15 GHz flux density (= 926 mJy). Shaded region: timespan of the Fermi peak (= (4± 1) X 1CT 6 ph cm~ 2 s -1 ), Dashed box: period of 
quiescent radio emission. 



gamma-ray peak without being observed. Ignoring the epoch 
beginnings and endings, which contain high flux values by 
construction, we find that the flux density of Cyg X-3 is larger 
than 250 mJy only 0.3% of the time during its quiescent state. 

5. SUMMARY & CONCLUSIONS 

We have presented observations of a 2010 May Cyg X-3 
flaring episode in the radio, soft and hard X-ray, and 7-ray 
bands. Our data show a sequence of three notable events: 
a particular softening and partial recovery of the X-ray emis- 
sion, a rapid ~1 Jy radio flare, and a 7-ray flare. Interpretation 
of the data is made more challenging by the sparse sampling 
of the radio data and the low cadence of the high-energy ob- 
servations as compared to the rapidity with which the radio 
emission can evolve. Nonetheless, the observations we do 
have and comparisons to historical data challenge the inter- 
pretation that the 7-ray emission is due to IC upscattering of 
the companion radiation field by high-energy electrons leav- 
ing the system in a relativistic jet, because there is no evidence 
for the subsequent radio emission that one would expect to see 
from these electrons. 

While the companion IC model of 7-ray emission from 
microquasars is clear and compelling, the lack of consistent 
radio/7-ray timing lags call into question its completeness. 
There could be more than one mechanism responsible for the 
7-ray emission of Cyg X-3, becoming more or less relevant in 
different circumstances, or the reprocessing of the 7-ray emis- 
sion by effects such as pair production could be more signifi- 
cant than commonly assumed. Detect ions of or limits to v ery 
high-energy 7-rays (>100 GeV; e.g. lAleksic et alj|2010l) or 
neutrinos from Cyg X-3 would aid in the understanding of the 
processes relevant to the emission in the Fermi band. 

The power of multiband monitoring of Cyg X-3 promises 



to increase significantly with the addition of 7-ray data to the 
set of available observations. As Cyg X-3 inevitably produces 
more flares, frequent observations at all wavelengths are im- 
portant to build a detailed understanding of the launching and 
propagation of its relativistic jets. Radio observations of <2 h 
cadence with near-continuous coverage would be ideal for es- 
tablishing an unambiguous relationship between radio and 7- 
ray flaring. Based on the phenomenology we observe, inten- 
sive radio observations should be triggered at the ends of very 
soft X-ray states without waiting for the detection of 7-ray 
flares. 
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